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1  Introduction 


Uninhabited  air  vehicles  (UAVs)  employed  for  reconnaissance  and  combat  missions  commonly  utilize  low- 
pressure  turbines  in  their  propulsion  systems.  Due  to  a  reduction  in  atmospheric  density  during  high-altitude 
cruise,  such  low-pressure  turbines  may  encounter  Reynolds  numbers,  based  upon  blade  axial  chord  and  inlet 
conditions,  below  25,000.  In  this  situation,  boundary  layers  along  a  large  extent  of  blade  surfaces  can  remain 
laminar,  even  in  the  presence  of  elevated  freestream  turbulence  levels.  The  laminar  boundary  layers  are  then 
particularly  susceptible  to  flow  separation  over  the  aft  portion  of  blade  suction  surfaces,  causing  blockage 
in  flow  passages  and  a  significant  reduction  in  turbine  efficiency.  These  circumstances  may  thus  result  in 
imposed  ceiling  limitations  for  prolonged  UAV  deployment. 

A  number  of  experimental  studies  [1,  2,  3,  4,  5,  6,  7,  8]  have  identified  several  factors  contributing 
to  performance  degradation  in  low-pressure  turbines,  which  in  addition  to  the  Reynolds  number  include 
freestream  turbulence,  pressure  gradient,  and  curvature.  Because  of  the  impact  upon  UAV  operations, 
a  number  of  experiments  have  focused  on  flow  control  as  a  means  of  reducing  efficiency  degradation  at 
altitude.  It  was  found  by  Lake  et  al. [9]  that  boundary-layer  separation  could  be  mitigated  by  modifying  the 
turbine  blade  geometry.  Dimples  were  recessed  into  blade  suction  surfaces  in  order  to  maintain  the  region 
of  attached  flow.  Although  this  passive  technique  increases  efficiency  at  low-Reynolds  numbers,  it  may  also 
increase  viscous  losses  at  higher  Reynolds  numbers  such  as  those  that  occur  during  takeoff,  where  unmodified 
blades  perform  satisfactorily. 

As  an  alternative,  Bons  et  al.  [10,  11,  12]  and  Sondergaard  et  al.  [13]  have  explored  active  flow  control  in 
the  form  of  both  steady  and  pulsed  vortex  generator  jets,  which  may  be  actuated  upon  demand.  Extensive 
measurements  verified  that  the  jets  dramatically  reduced  separation,  resulting  in  decreased  losses.  Plasma 
control  has  also  been  demonstrated  by  Huang  et  al.  [14]  to  restrict  separation.  More  recent  experiments  by 
Sondergaard  et  al.,  [15]  have  investigated  the  feasibility  of  increasing  the  blade  spacing  at  constant  chord 
for  a  linear  low-pressure  turbine  cascade,  thereby  raising  the  per  blade  loading.  For  practical  applications, 
a  higher  loading  can  reduce  the  turbine  part  count  and  stage  weight.  Increased  blade  spacing  however, 
is  accompanied  by  more  extensive  boundary-layer  separation  on  the  suction  surface  of  each  blade  due  to 
uncovered  turning,  resulting  in  a  further  reduction  of  efficiency  and  additional  wake  losses.  Vortex  generator 
jets  were  then  employed  to  diminish  these  losses  by  maintaining  attached  boundary-layer  flow  over  a  greater 
distance  along  the  blade  surface. 

Numerical  computations  of  low-pressure  turbines  have  traditionally  been  predicated  on  solution  of  the 
Reynolds-averaged  Navier-Stokes  (RANS)  equations.  [16,  17,  18,  19,  20]  Although  these  efforts  have  been 
moderately  successful  in  predicting  some  of  the  overall  features  of  the  turbine  flowfields,  limitations  and 
deficiencies  of  RANS  turbulence  models  have  precluded  an  accurate  description  of  low-Reynolds  separation 
and  the  transition  to  turbulence.  Moreover,  it  is  unlikely  that  the  RANS  approach  will  be  adequate  for 
simulating  active  flow  control  applications  in  the  transitional  regime. 

Due  to  advances  in  the  speed  and  storage  capacity  of  high-performance  computing  systems,  direct  numer¬ 
ical  simulation  (DNS)  and  large-eddy  simulation  (LES)  have  emerged  as  viable  means  for  the  investigation 
of  flows  through  low-pressure  turbines.  [21,  22,  23,  24,  25,  26,  27,  28]  DNS  and  LES  supported  by  high-order 
numerical  methods  improve  predictive  capability  where  RANS  models  fail,  and  allow  for  the  simulation  of 
flow  control  strategies.  Several  previous  computations  [23,  24,  27,  28]  have  been  related  to  the  experiments 
of  Bons  et  al.,  [10,  11,  12]  Sondergaard  et  al.,[15,  13]  and  Rivir  et  al. [29]  In  particular,  the  calculations 
of  Rizzetta  and  Visbal  [27]  considered  the  highly-loaded  low-pressure  turbine  configuration  investigated  by 
Sondergaard  et  al.  [15]  Simulations  were  carried  out  for  both  an  uncontrolled  baseline  case,  and  a  case 
with  active  flow  control  in  the  form  of  pulsed  vortex  generating  jets.  Although  the  aft-blade  separation  was 
mitigated  and  wake  total  pressure  losses  were  decreased  by  the  use  of  pulsed  injection,  the  jet  momentum 
coefficient  employed  in  the  computation  was  appreciably  lower  that  that  of  the  experiment.  Consequently, 
the  amount  of  control  attained  by  the  simulation  was  less  substantial  than  that  observed  by  Sondergaard 
et  al.  [15] 

Control  of  the  highly-loaded  low-pressure  turbine  considered  by  Sondergaard  et  al.  [15]  consisted  only 
of  steady  vortex  generator  jets.  And  while  prior  investigations  verified  that  separation  could  be  similarly 
managed  by  pulsed  injection,  [11,  12]  it  is  unclear  how  the  detailed  behavior  of  the  individual  flow  control 
arrangements  differs.  The  purpose  of  the  present  work  was  to  conduct  numerical  computations  of  the  turbine 
configuration  which  approximate  the  experiment  of  Sondergaard  et  al.  [15]  Steady  vortex  generator  jet  flow 


1 


control  matching  the  experimental  jet  momentum  coefficient  was  simulated  at  a  chord  inlet  Reynolds  number 
of  approximately  25,000.  In  addition,  a  calculation  was  also  performed  for  pulsed  jet  injection,  with  a  duty 
cycle  that  was  one  half  that  of  the  steady  case.  This  allows  for  a  direct  comparison  between  the  respective 
flow  control  situations,  as  well  as  with  the  baseline  case,  so  that  characteristics  of  the  flowfields  may  be 
elucidated. 
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2  The  Governing  Equations 


The  governing  equations  are  taken  as  the  unsteady  three-dimensional  compressible  unfiltered  Navier-Stokes 
equations.  After  introducing  a  curvilinear  coordinate  transformation  to  a  body- fitted  system,  the  equations 
are  cast  in  the  following  nondimensional  conservative  form 
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=  0. 


(1) 


Here  t  is  the  time,  £,  77,  £  the  computational  coordinates,  Q  the  vector  of  dependent  variables,  F,  G,  H  the 
inviscid  flux  vectors,  and  Fv ,  Gv,  Hv  the  viscous  flux  vectors.  The  vector  of  dependent  variables  is  given  as 
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and  the  vector  fluxes  by 
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In  the  preceding  expressions,  u,v,w  are  the  Cartesian  velocity  components,  p  the  density,  p  the  pressure, 
and  T  the  temperature.  All  length  scales  have  been  nondimensionalized  by  the  turbine  blade  axial  chord 
c,  and  dependent  variables  have  been  normalized  by  their  reference  values  except  for  p  which  has  been 
nondimensionalized  by  Poo^^,-  Components  of  the  stress  tensor  and  heat  flux  vector  are  expressed  as 
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The  Sutherland  law  for  the  molecular  viscosity  coefficient  p  and  the  perfect  gas  relationship 
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were  also  employed,  and  Stokes’  hypothesis  for  the  bulk  viscosity  coefficient  has  been  invoked. 
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3  The  Numerical  Method 


Time-accurate  solutions  to  Eq.  (1)  are  obtained  numerically  by  the  implicit  approximately-factored  finite- 
difference  algorithm  of  Beam  and  Warming  [30]  employing  Newton- like  subiterations,  [31]  which  has  evolved 
as  an  efficient  tool  for  generating  solutions  to  a  wide  variety  of  complex  fluid  flow  problems,  and  may  be 
written  as  follows 
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In  this  expression,  which  is  employed  to  advance  the  solution  in  time,  Qp+1  is  the  p  +  1  approximation 
to  Q  at  the  n  +  1  time  level  Qn+1,  and  A  Q  =  Qp+1  —  Qp.  For  p  =  1,  Qp  =  Qn.  Second-order-accurate 
backward-implicit  time  differencing  was  used  to  obtain  temporal  derivatives. 

The  implicit  segment  of  the  algorithm  incorporates  second-order-accurate  centered  differencing  for  all 
spatial  derivatives,  and  utilizes  nonlinear  artificial  dissipation [32]  to  augment  stability.  Efficiency  is  enhanced 
by  solving  this  implicit  portion  of  the  factorized  equations  in  diagonalized  form.  [33]  Temporal  accuracy,  which 
can  be  degraded  by  use  of  the  diagonal  form,  is  maintained  by  utilizing  subiterations  within  a  time  step. 
This  technique  has  been  commonly  invoked  in  order  to  reduce  errors  due  to  factorization,  linearization, 
diagonalization,  and  explicit  application  of  boundary  conditions.  It  is  useful  for  achieving  temporal  accuracy 
on  overset  zonal  mesh  systems,  and  for  a  domain  decomposition  implementation  on  parallel  computing 
platforms.  Any  deterioration  of  the  solution  caused  by  use  of  artificial  dissipation  and  by  lower-order  spatial 
resolution  of  implicit  operators  is  also  reduced  by  the  procedure.  Three  subiterations  per  time  step  have 
been  applied  to  preserve  second-order  temporal  accuracy. 

The  compact  difference  scheme  employed  on  the  right-hand  side  of  Eq.  10  is  based  upon  the  pentadiagonal 
system  of  Lele,[34]  and  is  capable  of  attaining  spectral-like  resolution.  This  is  achieved  through  the  use  of  a 
centered  implicit  difference  operator  with  a  compact  stencil,  thereby  reducing  the  associated  discretization 
error.  For  the  present  computations,  a  fourth-order  tridiagonal  subset  of  Lele’s  system  is  utilized,  which  is 
illustrated  here  in  one  spatial  dimension  as 


a 


(11) 


with  a  =  1/4  and  a  =  3/2.  The  scheme  has  been  adapted  by  Visbal  and  Gaitonde[35]  as  an  implicit  iterative 
time-marching  technique,  applicable  for  unsteady  vortical  flows.  It  is  used  in  conjunction  with  a  sixth-order 
low-pass  Pade-type  non-dispersive  spatial  filter  developed  by  Gaitonde  et  al., [36]  which  has  been  shown  to 
be  superior  to  the  use  of  explicitly  added  artificial  dissipation  for  maintaining  both  stability  and  accuracy 
on  stretched  curvilinear  meshes. [35]  The  filter  is  applied  to  the  solution  vector  sequentially  in  each  of  the 
three  computational  directions  following  each  subiteration,  and  is  implemented  in  one  dimension  as 
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where  Q  designates  the  filtered  value  of  Q.  Equation  12  represents  a  one-parameter  family  of  sixth-order 
filters,  where  numerical  values  for  a /  and  the  a„’s  may  be  found  in  Ref.  Gaitonde: AFRL3060 

The  aforementioned  features  of  the  numerical  algorithm  are  embodied  in  a  parallel  version  of  the  time- 
accurate  three-dimensional  computer  code  FDL3DI,[37]  which  has  proven  to  be  reliable  for  steady  and 
unsteady  fluid  flow  problems,  including  the  simulation  of  flows  over  delta  wings  with  leading-edge  vortices,  [38, 
39,  40]  vortex  breakdown,  [39,  40]  direct  numerical  simulation  of  transitional  wall  jets [41]  and  synthetic  jet 
actuators,  [42]  and  DNS  and  LES  of  subsonic[24,  43]  and  supersonic  flowfields.[44,  45,  46] 
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4  Details  of  the  Computations 

The  experiment  Sondergaard  et  al.  [13]  consisted  of  a  linear  cascade  of  turbine  blades  having  a  span  of  0.88  m 
and  an  axial  chord  of  0.089  m.  The  axial  chord  to  spacing  ratio  (solidity)  was  0.75,  the  inlet  flow  angle 
a*  =  35.0  deg,  and  the  design  exit  flow  angle  a0  =  60.0  deg.  Shown  in  Fig.  1  is  a  schematic  representation 
of  the  turbine  blade  shape,  given  by  the  Pratt  &  Whitney  “PakB”  research  design,  which  is  a  Mach  number 
scaled  version  of  geometries  typically  used  in  highly-loaded  low-pressure  turbines.  [10,  11,  12,  15,  13]  These 
experiments  investigated  the  use  of  vortex  generator  jets  to  control  the  flow  by  mitigating  the  effects  of 
separation  on  the  blade  suction  surface.  The  jets  were  created  by  blowing  air  through  holes  which  had  been 
drilled  in  the  blade  surface  at  a  pitch  angle  of  30  deg  and  a  skew  angle  of  90  deg.  Here,  the  pitch  is  defined 
as  the  angle  the  jet  makes  with  the  local  surface,  and  the  skew  is  the  angle  of  the  projection  of  the  jet  on 
the  surface,  relative  to  the  “local  freestream”  direction.  [12]  The  size  of  the  drill  used  to  develop  the  holes  is 
commonly  referred  to  as  the  jet  diameter  d,  which  was  0.001  m.  Because  of  the  orientation  however,  the  jet 
exit  geometric  shape  is  elliptic  as  seen  in  Fig.  2,  and  the  jet  exit  velocity  vector  has  components  only  in  the 
blade-normal  (n)  and  spanwise  (z)  directions. 

4.1  Computational  Meshes 

To  conserve  computational  resources,  only  a  single  turbine  blade  passage  is  considered,  and  periodic  condi¬ 
tions  are  enforced  in  the  vertical  direction  (y)  to  represent  a  single  turbine  stage  flowfield.  This  situation  is 
somewhat  different  than  the  experimental  arrangement,  which  consisted  of  a  small  finite  number  of  blades 
in  a  linear  cascade.  The  computational  domain  surrounding  the  blade  was  described  by  a  body-fitted  mesh 
system,  whose  origin  was  located  at  the  inboard  leading  edge  of  the  blade  (see  Fig.  1).  The  mesh  employs  an 
O-grid  topology,  and  was  elliptically  generated  using  automated  software.  [47]  Figure  3  indicates  the  baseline 
grid  about  the  turbine  blade,  which  was  comprised  of  348  points  in  the  circumferential  direction  (J),  189 
points  in  the  blade-normal  direction  (J),  and  101  points  in  the  spanwise  direction  (K).  Minimum  spacing 
in  the  J-direction  occurs  at  the  blade  surface.  Mesh  points  for  1  <  I  <  5  and  Jmax_ 4  <  I  <  /max  are 
coincident  in  an  overlap  region  at  the  blade  leading  edge  so  that  periodic  conditions  in  the  circumferential 
direction  may  be  enforced  in  order  to  complete  the  O-grid  construct.  To  accommodate  periodic  conditions 
in  the  vertical  direction,  the  streamwise  (a:)  point  locations  and  vertical  spacing  along  /lu  —  J2tl  is  identical 
to  that  along  In  1 21  for  t/max— 4  —  J  —  /max- 

To  properly  capture  the  correct  fluid  physics  for  the  flow  control  cases,  the  baseline  grid  was  modified  in 
order  to  enhance  resolution  in  the  jet  and  jet  wake  regions.  This  was  done  by  removing  points  in  the  near 
wall  grid,  indicated  by  the  blanked  area  of  Fig.  3,  and  replacing  them  with  an  embedded  refined  mesh  as 
seen  in  Fig.  4.  This  technique  is  identical  to  that  which  was  successfully  employed  by  Rizzetta  and  Visbal[27] 
in  a  similar  computation.  The  size  of  the  refined  mesh  region  was  (313  x  185  x  101)  grid  points  in  (/,  J,  K) 
respectively.  A  series  of  overset  grids  were  used  to  describe  the  elliptic  jet  geometry,  and  zones  upstream, 
downstream,  inboard,  and  outboard  of  the  jet.  These  grids  connect  to  each  other  and  with  the  refined-mesh 
region  as  shown  in  Fig.  5.  Also  seen  in  the  figure  is  the  jet  nozzle  interior  to  the  blade,  which  was  accounted 
for  in  the  simulations,  and  extended  a  distance  of  four  jet  diameters  below  the  surface. 

In  order  to  facilitate  application  of  inflow  and  outflow  conditions  to  the  turbine  blade  domain,  overset 
grids  were  utilized  upstream  and  downstream  of  the  blade  region.  These  are  found  in  Fig.  6,  and  consisted 
of  (46  x  61  x  101)  and  (31  x  75  x  101)  mesh  points  in  he  streamwise,  vertical,  and  spanwise  directions  for  the 
upstream  and  downstream  domains  respectively.  Half  the  grid  points  in  the  vertical  direction  are  seen  in  the 
figure.  Although  the  experiment  of  Sondergaard  et  al.  [15]  consisted  of  a  row  of  jets  uniformly  distributed 
along  the  span  of  a  blade,  only  a  single  jet  was  considered  in  the  simulations.  The  z-extent  of  the  domain 
was  taken  as  the  inter-jet  spacing  s  =  0.112,  which  is  equal  to  10  jet  diameters  (as  in  the  experiment)  and 
periodic  conditions  were  invoked  at  the  spanwise  boundaries.  A  non-uniform  spanwise  mesh  distribution  was 
utilized,  with  grid  point  clustering  near  the  jet  exit.  At  all  overset  domain  boundaries,  a  minimum  5-point 
overlap  was  constructed.  The  total  number  of  grid  points  in  all  domains  was  approximately  13.5  x  106. 
Computational  mesh  spacings  given  in  wall  units  for  both  simulations  are  presented  in  Table  1.  Results  of 
the  grid  refinement  study  performed  for  the  baseline  uncontrolled  case  of  Ref.  Rizzetta:2004-2204,  indicates 
that  adequate  resolution  should  be  achieved  for  the  present  computations. 
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Table  1:  Computational  mesh  spacings  in  wall  units. 


Case 

AC1" 

A4in 

A /+ 

^''max 

A4in 

A4in 

steady 

0.02 

1.63 

50.07 

0.44 

3.09 

10.98 

pulsed 

0.02 

1.74 

53.39 

0.46 

3.23 

11.47 

4.2  Temporal  Considerations 

Calculations  for  the  steady  jet  case  were  obtained  using  a  time  step  of  At  =  1.50  x  10-4,  where  t  is  nondi- 
mensionalized  by  reference  quantities.  Based  upon  inflow  conditions,  this  corresponds  to  a  time  increment  of 
approximately  1.18  x  10~4  (see  Table  1  for  the  time  step  in  wall  units).  For  the  pulsed  case,  these  parameters 
were  At  =  1.47  x  10-4  and  At,  =  1.20  x  10-4  respectively.  A  slightly  smaller  time  step  was  employed  with 
unsteady  injection  so  that  the  pulsing  cycle  could  be  described  by  an  integer  number  of  increments  (1300 
time  steps/cycle).  Flowfields  for  each  case  were  initialized  from  the  previously  obtained  baseline  (uncon¬ 
trolled)  solution,  and  processed  for  20,000  time  steps  in  order  to  remove  transients.  Final  results  were  then 
evolved  for  a  subsequent  124,000  steps,  during  which  statistical  information  was  collected.  Statistics  were 
monitored  to  assure  that  a  converged  sample  was  achieved.  This  duration  represents  over  14.5  characteristic 
time  units,  based  upon  the  inflow  velocity,  for  both  simulations. 

Experimental  investigations[10,  11,  12,  15,  13]  have  used  several  parameters  to  describe  characteristics 
of  vortex  generating  jets.  The  jet  blowing  ratio  B  is  defined  as 

B  =  —  (13) 

Qi 

where  qj  is  the  jet  exit  velocity  magnitude,  and  is  the  “local  freestream  velocity”  (see  Ref.  Bons:JTM2001 
for  the  definition  of  g;).  In  the  present  computations,  the  maximum  value  of  the  jet  blowing  ratio  _Bmax  was 
specified  as  4.5.  This  is  identical  to  the  experiment  of  Sondergaard  et  al. [15]  for  the  steady  case.  As  a  point 
of  reference,  typical  film  cooling  holes  for  turbine  applications  are  designed  to  operate  at  an  approximate 
value  of  F>max  =  2.0. [12]  For  the  present  configuration  considered  here,  experimental  data  is  available  only 
for  steady  blowing  flow  control.  In  previous  investigations  with  a  smaller  inter-blade  spacing  however, 
measurements  were  also  taken  using  pulsed  mass  injection.  The  unsteady  jet  exit  velocity  sequence  in  those 
experiments  consisted  of  a  period  of  essentially  constant  blowing,  followed  by  an  interval  of  no  blowing. 
The  duty  cycle  represents  that  portion  of  the  total  duration  over  which  the  jet  is  active,  and  was  50% 
for  the  present  case.  A  nondimensional  frequency  T  as  formulated  in  Refs.  Bons:99-0367,  Bons:JTM2001, 
Bons:JTM2002,  Sondergaard:JPP2002  was  used  to  quantify  the  forcing  period  where 

T  =  OAlcf/u^q,.  (14) 

The  value  of  T  =  3.1  was  used  in  the  present  computations,  which  was  the  same  as  that  of  some  experiments[10, 
11,  12].  The  jet  blowing  coefficient,  defined  as 

c^  =  V  [(tt/4 )d2£Lx]  (15) 

is  used  to  quantify  the  momentum  flux  transfered  to  the  turbine  blade  passage  by  active  control.  For  the 
steady  blowing  case  =  2.0  x  10-3,  and  is  one  half  this  value  for  pulsed  injection  owing  to  the  50%  duty 
cycle. 

A  representative  experimental[12]  time  history  of  one  blowing  cycle  obtained  by  hot  wire  measurement 
is  shown  in  Fig.  7.  For  the  pulsed-injection  computation,  a  series  of  cubic  and  linear  functions  were  used  to 
approximate  the  experimental  result,  as  seen  in  the  figure  where  the  50%  duty  cycle  is  apparent.  As  noted 
earlier,  a  single  cycle  was  described  by  1300  time  steps,  and  the  total  computation  consisted  of  96  cycles. 

4.3  Boundary  Conditions 

Inflow  and  outflow  conditions  for  the  complete  turbine  blade  domain  were  obtained  in  a  manner  consistent 
with  subsonic  internal  flows,  which  have  commonly  been  employed  in  RANS  applications. [48]  Along  the 


6 


upstream  boundary,  the  total  pressure,  total  temperature,  and  inlet  flow  angle  were  specified,  and  the 
velocity  magnitude  was  obtained  from  the  interior  solution.  Downstream,  the  exit  static  pressure  (poo) 
was  fixed,  and  other  flow  variables  were  extrapolated  from  within  the  domain.  Because  the  inflow  velocity 
develops  as  part  of  the  solution,  the  Reynolds  number  based  upon  the  inflow  conditions  was  not  known  a 
priori.  A  reference  Mach  number  (M0 c)  of  0.1  was  selected  for  all  computations.  For  the  baseline  solution, 
the  reference  Reynolds  number  (Re oo)  was  adjusted  to  match  the  desired  inflow  condition,  Re.  This  same 
Reynolds  number  was  then  employed  in  the  current  simulations.  When  flow  control  was  applied  however, 
the  inflow  velocity  magnitude  increased  due  to  a  reduction  of  blockage  in  the  blade  passage.  Although 
this  situation  was  anticipated,  the  reference  Reynolds  number  was  not  altered.  Thus,  the  inflow  Reynolds 
numbers  are  somewhat  higher  than  the  nominal  value  of  25,000.  For  the  time-mean  flowfields,  Re  =  29, 100 
for  steady  blowing  and  Re  =  29,960  with  pulsed  mass  injection.  This  approach  for  the  inflow  and  outflow 
boundaries  was  previously  employed  by  Rizzetta  and  Visbal[24,  27]  in  similar  computations. 

Periodic  conditions  were  applied  along  the  upper  and  lower  portions  of  the  turbine  blade  as  indicated 
previously  in  Fig.  3.  Periodic  conditions  were  also  applied  along  upper  and  lower  horizontal  boundaries  of 
the  upstream  and  downstream  domains  seen  in  Fig.  6.  The  downstream  domain  was  intensionally  severely 
stretched  in  the  streamwise  direction  to  prevent  spurious  reflections  from  the  outflow  boundary.  This  tech¬ 
nique  transfers  information  to  high  spatial  wave  numbers,  and  then  dissipates  it  by  the  low-pass  filter.  [49] 
Flow  variables  in  all  regions  of  overlapped  meshes  were  obtained  from  explicit  sixth-order  accurate  La- 
grangian  interpolation  formulae.  This  includes  the  upstream  and  downstream  domains,  the  turbine  blade 
grid,  the  refined-mesh  region,  all  zones  surrounding  the  jet,  and  the  interior  jet  nozzle.  The  interpolation 
approach  for  high-order  numerical  solutions  has  been  successfully  applied  by  Sherer[50]  for  the  simulation  of 
fluid  dynamic  and  acoustic  problems.  As  stated  previously,  periodic  conditions  were  invoked  in  the  spanwise 
direction. 

On  solid  surfaces,  including  the  blade  boundary  and  walls  of  the  jet  nozzle,  the  no  slip  condition  was 
enforced  along  with  an  isothermal  wall,  and  a  fourth-order  accurate  representation  of  zero  normal  pressure 
gradient.  At  the  inflow  to  the  jet  nozzle,  a  velocity  profile  was  specified,  the  jet  was  assumed  to  be  isothermal, 
and  the  pressure  was  obtained  from  the  inviscicl  normal  momentum  equation.  For  both  the  steady  and  pulsed 
cases,  the  velocity  profile  was  taken  to  be  parabolic.  In  the  steady  case,  the  maximum  jet  blowing  ratio  was 
fixed  at  Bm ax  =  4.5,  corresponding  to  the  experiment,  while  for  pulsed  injection  this  ratio  varied  in  time 
according  to  the  form  shown  in  Fig.  7.  It  is  noted  for  both  cases  that  the  maximum  jet  velocity  is  over  seven 
time  greater  than  the  magnitude  of  the  blade  inflow  velocity.  The  Reynolds  number  of  the  jet  based  upon 
the  maximum  blowing  ratio  and  the  diameter  d  was  2370. 

4.4  Domain  Decomposition 

For  parallel  processing,  the  previously  described  computational  domains  were  decomposed  into  a  series  of 
subzones,  which  were  then  distributed  on  individual  processors  of  a  massively  parallel  computing  platform 
(IBM  SP4,  SGI  3900).  Decompositions  were  constructed  to  provide  an  approximately  equal  number  of  grid 
points  in  every  subzone,  thereby  balancing  the  computational  work  load  among  the  processors.  Faces  at  the 
boundaries  of  each  subzone  block  were  overset  into  adjacent  domains,  such  that  an  overlap  of  five  planes  was 
established.  Although  this  incurred  an  overhead  due  to  redundant  computation,  it  maintained  the  formal 
high-order  accuracy  of  both  the  numerical  differencing  and  filtering  schemes.  Because  a  vast  majority  of 
the  overlapping  mesh  points  of  the  respective  decompositions  were  coincident,  no  further  interpolation  was 
required.  Automated  software[51]  was  used  to  identify  donor  and  recipient  grid  points  in  the  overlapping 
domains.  Inter-node  communication  among  the  processors  was  established  through  standard  message-passing 
interface  (MPI)  library  routines, [52]  which  were  used  to  transfer  information  between  the  various  subzones 
of  the  flowfield  at  domain  boundaries.  A  total  of  202  processors  were  employed  for  all  of  the  computations 
reported  here. 
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5  Results 

5.1  The  Time-Mean  Flowfields 


Time-mean  surface  pressure  coefficient  distributions  for  the  flow  control  and  baseline  [27]  computations  are 
found  in  Fig.  8.  For  the  flow  control  cases,  these  distributions  were  obtained  at  the  spanwise  location  of 
the  periodic  boundary  between  two  jets.  As  indicated  in  the  figure,  the  jet  was  positioned  at  x  =  0.37, 
which  coincides  with  the  separation  point  of  the  baseline  solution.  The  large  plateau  region  in  the  baseline 
distribution  is  characteristic  of  massively  separated  flow.  Because  of  reduced  blockage  and  increased  inflow 
velocity,  the  effect  of  flow  control  is  to  decrease  the  pressure  on  the  upstream  portion  of  the  suction  surface, 
while  increasing  it  downstream,  relative  to  the  baseline  case.  Separation  occurs  at  x  =  0.56  for  the  steady 
blowing  result,  and  at  x  =  0.58  for  pulsed  injection.  Only  minor  differences,  most  noticeable  near  the  jet, 
are  observed  between  the  two  flow  control  solutions.  Also  shown  in  the  figure,  is  a  coarse-mesh  distribution 
for  the  baseline  case,  which  indicates  little  sensitivity  to  grid  resolution,  except  near  the  trailing-edge  region. 
More  extensive  results  of  the  grid  resolution  study  for  the  baseline  case  are  found  in  Ref.  Rizzetta:2004-2204. 

Profiles  of  the  time-mean  velocity  magnitude  are  presented  in  Figs.  9  and  10.  These  profiles  were  extracted 
along  lines  (n)  normal  to  the  blade  surface  at  each  a:-station,  on  the  periodic  spanwise  boundary  between 
jets.  The  data  of  Sondergaard  et  al.  [15]  for  steady  blowing  also  appears  in  the  figures.  However,  only  a 
qualitative  comparison  between  the  experimental  and  numerical  results  may  be  expected  due  to  appreciable 
differences  in  the  respective  flow  situations.  Among  these  are  the  fact  that  the  experiments  were  conducted  in 
a  linear  cascade,  which  had  a  limited  number  of  blades  because  of  the  large  inter-blade  spacing.  Tail  boards 
were  attached  at  the  trailing  edge  of  the  first  and  last  blades  in  the  cascade,  and  extended  downstream  a 
considerable  distance  into  the  wake,  thus  creating  covered  turning  for  a  substantial  portion  of  the  flowfield. 
Vortex  generating  jets  were  fabricated  only  on  the  central  blade  of  the  cascade.  This  arrangement  is  quite 
unlike  the  simulations,  where  periodic  conditions  about  the  blade  are  enforced.  The  computations  were 
carried  out  at  a  nominal  Reynolds  number  of  25,000  which  is  representative  of  the  practical  situation  for 
UAVs,  and  may  be  compared  with  the  previously  obtained  baseline  case.  Experimental  data  was  collected 
for  Re  =  50,  000. 

Profiles  in  Figs.  9  and  10  have  been  normalized  by  the  inflow  velocity  magnitude.  The  major  effect  of 
flow  control  is  to  increase  fullness  of  the  profiles  through  momentum  transfer  from  the  jet  and  mixing  with 
the  boundary  layer.  Only  minor  variations  are  observed  between  the  steady  and  pulsed  injection  cases. 
As  expected,  profiles  from  the  simulations  are  thicker  than  those  of  the  experiment  due  to  aforementioned 
differences  in  the  Reynolds  numbers  and  configurations.  Increased  fullness  of  the  control  profiles  from  the 
baseline  solution  is  evident  in  the  figures,  even  at  locations  where  separation  is  present  (a:  >  0.58). 

Spanwise  turbulent  kinetic  energy  wave-number  spectra  are  provided  in  Fig.  11.  These  spectra  were 
generated  along  lines  in  the  ^-direction  at  a  distance  of  n  =  0.03  from  the  blade  surface.  This  location 
(n  =  0.03)  is  approximately  equal  to  one  half  of  the  boundary-layer  thickness  of  the  time-mean  velocity 
profile  upstream  of  separation.  At  the  two  most  upstream  stations  (x  =  0.50,0.70),  E^z  is  higher  for  the 
flow  control  cases  than  that  of  the  baseline  due  to  energy  being  added  to  the  flow.  In  the  downstream 
region  (x  =  0.90),  flow  control  has  mitigated  separation  and  the  associated  breakdown  into  a  more  chaotic 
situation,  so  that  the  turbulent  kinetic  energy  of  the  baseline  case  is  higher.  Because  of  the  low  Reynolds 
number,  only  a  small  portion  of  the  spectrum  at  low  wave  number  lies  in  the  inertial  range. 

Time-mean  contours  of  the  streamwise  component  of  velocity  u,  are  seen  in  Fig.  12.  The  contours  for 
the  flow  control  results  were  taken  at  the  plane  of  the  periodic  boundary  between  jets,  while  those  of  the 
baseline  case  have  been  spanwise  averaged.  An  increased  region  of  attached  flow  due  to  flow  control  is 
apparent,  resulting  in  a  decrease  of  the  wake  thickness.  The  effect  of  the  decrease  is  quantified  by  the 
integrated  wake  total  pressure  loss  coefficient  Cw,  defined  as[15] 

C"  =  s/(S)^  (16) 

where  the  integration  takes  place  across  the  wake  profile  between  the  periodic  boundaries  in  the  vertical 
direction.  It  was  determined  that  Cw  =  1.10 [27]  for  the  baseline  solution,  Cw  =  0.52  for  steady  blowing, 
and  Cw  =  0.48  for  pulsed  injection.  Thus  even  with  a  50%  duty  cycle,  pulsed  injection  was  found  to 
be  as  effective  as  steady  blowing  for  the  reduction  of  wake  losses.  This  result  is  consistent  with  previous 


experimental  studies.  [11,  12] 

Corresponding  time-mean  contours  of  the  spanwise  component  of  vorticity  are  displayed  in  Figure  13. 
Maintaining  attached  flow  and  decreasing  the  extent  of  the  wake  are  again  exhibited  by  the  flow  control 
results.  Fig.  14  represents  these  same  contours  evaluated  on  the  blade  surface  and  viewed  from  above.  Apart 
from  the  area  immediately  surrounding  the  jet,  the  flow  control  cases  are  remarkably  similar. 

5.2  Unsteady  Features  of  the  Flowfields 

Profiles  of  the  root-mean-square  fluctuating  velocity  magnitude  at  several  steamwise  locations  are  found  in 
Fig.  15.  At  the  upstream  stations  (x  =  0.50,0,70),  the  flow  control  cases  differ  somewhat  from  each  other, 
but  both  produce  peak  values  less  than  that  of  the  baseline  result.  Fluctuations  in  the  control  cases  are 
related  to  unsteadiness  of  the  jet,  while  those  of  the  baseline  flow  evolve  from  the  shear  layer.  As  a  result, 
the  location  of  the  maximum  fluctuation  occurs  further  from  the  blade  surface  in  the  baseline  solution.  At 
the  downstream  location  ( x  =  0.90),  a  more  massive  separated  region  in  the  baseline  flow  gives  rise  to  a 
uniformly  high  fluctuating  velocity  in  the  normal  direction,  which  is  more  localized  in  the  flow  control  cases. 

Turbulent  kinetic  energy  frequency  spectra  are  shown  in  Fig.  16.  The  data  used  to  generate  these  spectra 
was  collected  at  n  =  0.03,  similar  to  that  of  Fig.  11.  Like  the  spanwise  wave  number  spectra  of  Fig.  11,  the 
frequency  spectra  for  the  flow  control  cases  is  substantially  higher  than  that  of  the  baseline  result  at  the 
upstream  stations  (x  =  0.50,0.70)  because  of  energy  addition.  At  x  =  0.50,  the  discrete  peaks  in  Eu  for  the 
pulsed  injection  case  correspond  to  harmonics  of  the  forcing  frequency.  The  occurrence  of  this  behavior  was 
described  in  detail  in  Ref.  Rizzetta:2004-2204.  As  noted  previously,  turbulent  energy  of  the  baseline  case 
exceeds  that  of  the  flow  control  results  near  the  blade  trailing  edge  (x  =  0.90). 

Instantaneous  contours  of  the  streamwise  velocity  component  u  are  presented  Fig.  17.  The  baseline  result 
was  extracted  at  the  midspan  location,  while  those  of  the  flow  control  cases  were  situated  between  control 
jets.  Similar  to  the  time-mean  contours  of  Fig.  12,  it  is  noted  that  flow  control  maintains  attached  flow 
and  decreases  the  vertical  extent  of  the  wake  relative  to  the  baseline  case.  Unsteady  structures  are  visible 
in  separated  flow  regions.  These  structures  are  more  easily  discerned  in  instantaneous  spanwise  vorticity 
contours  appearing  in  Fig.  18.  It  is  seen  in  the  flow  control  cases  that  vorticity  is  being  generated  in  the 
boundary  layer  in  close  proximity  to  the  blade  surface.  When  the  boundary  layers  separate  in  these  situations, 
it  is  much  less  dramatic  than  in  the  uncontrolled  case.  For  the  baseline  solution,  the  extensive  unsteady 
separated  flow  region  has  a  richer  content  of  small  scale  structures  due  to  breakdown  and  a  transition  to  a 
more  complex  situation.  A  three-dimensional  representation  of  the  flow  is  depicted  by  iso-surfaces  of  vorticity 
magnitude  in  the  trailing-edge  region  seen  in  Fig.  19.  The  value  of  iso-surfaces  correspond  approximately 
to  that  at  the  edge  of  the  shear  layer  upstream  of  separation.  Both  the  vertical  and  spanwise  extent  of  the 
turbulent  structures  are  visible  for  each  case. 

Instantaneous  spanwise  vorticity  contours  on  the  blade  surface  viewed  from  above  are  displayed  in  Fig.  20. 
Although  differences  in  the  control  case  exist  in  the  near-jet  region,  the  downstream  flowfields  are  quite 
similar.  Iso-surfaces  of  vorticity  magnitude  in  the  near-jet  region,  as  viewed  downstream  of  the  jet,  are 
provided  in  Fig.  21.  Unsteadiness  can  be  observed  in  the  jet  surface,  just  downstream  of  the  nozzle  exit,  for 
the  steady  blowing  case.  It  is  noted  that  the  wake  of  the  jet  extends  across  the  spanwise  domain.  The  jet 
surface  for  the  pulsed  case  is  less  coherent. 

Instantaneous  planar  contours  of  the  streamwise  component  of  vorticity  are  found  in  Figs.  22-24.  Lo¬ 
cations  of  these  contours  lie  at  distances  of  1,  5,  and  10  jet  diameters  downstream  from  the  center  of  the 
jet.  The  height  of  the  attached  boundary  layer  (ft)  upstream  of  the  jet  is  indicated  in  the  figures.  For  the 
steady  injection  case  of  Fig.  22,  the  jet  penetrates  to  a  height  greater  than  the  boundary-layer  thickness, 
and  mixing  persists  downstream  beyond  x  —  Xj  =  10 d.  In  Fig.  23,  contours  are  shown  for  the  pulsed  case,  at 
a  time  corresponding  to  the  midpoint  of  the  pulsing  cycle  (t/tp  =  0.5).  At  this  instance,  the  jet  exit  velocity 
is  at  its  maximum  and  has  been  active  for  one  half  the  cycle  (see  Fig.  7).  Here,  the  contours  look  similar  to 
those  of  the  steady  injection  case,  except  at  the  most  downstream  station.  Contours  for  the  pulsed  case  at 
the  end  of  the  injection  cycle  appear  in  Fig.  24.  For  this  situation  (t/tp  =  1.0),  no  injection  has  occurred  for 
one  half  of  the  cycle,  and  effects  of  the  jet  are  only  visible  at  the  most  downstream  station. 

It  is  perhaps  surprising  that  there  are  so  few  differences  between  the  steady  blowing  and  pulsed  injection 
cases.  The  reason  for  this  is  that  even  in  the  steady  case,  the  flow  immediately  downstream  of  the  jet  is 
inherently  unsteady.  Because  of  the  orientation  of  the  jet  transverse  to  the  streamwise  flow  (blade-normal 
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and  spanwise),  it  is  fundamentally  unstable  and  rapidly  evolves  to  a  non-coherent  form.  This  process  serves 
to  generate  mixing,  much  like  the  pulsed  injection  case.  The  pulsed  case  however,  provides  a  much  more 
efficient  means  of  distributing  momentum  to  the  flow,  since  it  produces  identical  benefits  with  one  half  the 
mass  flow  rate. 
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6  Summary  and  Conclusion 

Direct  numerical  simulation  was  used  to  describe  the  transitional  subsonic  flow  through  a  highly-loaded  low- 
pressure  turbine  at  a  nominal  Reynolds  number  of  25,000,  with  both  steady  blowing  and  pulsed  mass  injection 
flow  control.  Computations  were  carried  out  using  a  high-order  numerical  method,  and  an  overset  grid 
methodology  for  describing  the  computational  flowfield.  A  high-order  interpolation  procedure  maintained 
accuracy  at  boundaries  of  discrete  mesh  systems  that  were  used  to  define  the  jet  region,  including  its  interior 
nozzle. 

For  both  the  steady  and  pulsed  cases,  the  fundamental  effect  of  flow  control  was  to  energize  the  blade 
boundary  layer  due  to  the  transfer  of  fluid  momentum  and  mixing.  This  helped  maintain  attached  flow  along 
the  blade  surface  for  a  distance  of  19-21%  greater  than  that  of  the  baseline  case.  As  a  result,  wake  total 
pressure  losses  were  decreased  by  53-56%. 

Comparisons  were  made  for  both  the  time-mean  and  fluctuating  components  of  the  respective  flowfields. 
Velocity  magnitude  profiles  qualitatively  agreed  with  experimental  measurements  that  were  taken  for  a  some¬ 
what  different  configuration  at  a  higher  Reynolds  number.  Few  differences  were  observed  in  characteristics 
of  the  flow  control  solutions.  Although  the  steady  blowing  case  employed  no  explicit  unsteady  forcing,  due 
to  fluid  instability  of  the  jet,  the  flowfield  was  inherently  unsteady  and  resulted  in  features  which  were  quite 
similar  to  the  pulsed  case.  Both  situations  achieved  the  desired  result  of  decreased  wake  losses,  but  the 
pulsed  case  required  less  mass  flow  to  attain  the  same  level  of  improvement  in  efficiency. 


11 


Figure  1:  Schematic  representation  of  the  turbine  blade  configuration. 
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Figure  3:  Turbine  blade  computational  mesh. 
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Figure  4:  Refined-mesh  region. 


Figure  5:  Vortex-generator  jet  region. 
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Figure  6:  Upstream  and  downstream  computational  meshes. 
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Figure  7:  Jet  blowing  ratio  time  history. 
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Figure  8:  Time-mean  surface  pressure  coefficient  distributions. 


Figure  9:  Time-mean  velocity  magnitude  profiles  at  upstream  stations. 
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Figure  10:  Time-mean  velocity  magnitude  profiles  at  downstream  stations. 
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Figure  11:  Time-mean  turbulent  kinetic  energy  spanwise  wave-number  spectra. 
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Figure  12:  Time-mean  planar  contours  of  u  velocity. 


23 


baseline 


steady 


pulsed 


Figure  13:  Time-mean  contours  of  the  spanwise  component  of  vorticity. 
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Figure  14:  Time-mean  contours  of  the  spanwise  component  of  vorticity  on  the  blade  surface 
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Figure  15:  Time-mean  fluctuating  velocity  magnitude  profiles. 
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Figure  16:  Turbulent  kinetic  energy  frequency  spectra. 
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Figure  17:  Instantaneous  planar  contours  of  u  velocity. 
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Figure  18:  Instantaneous  planar  contours  of  the  spanwise  component  of  vorticity. 
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Figure  19:  Instantaneous  iso-surfaces  of  vorticity  magnitude  in  the  trailing-edge  region. 
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Figure  20:  Instantaneous  planar  contours  of  the  spanwise  component  of  vorticity  on  the  blade  surface. 
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Figure  21:  Instantaneous  iso-surfaces  of  vorticity  magnitude  in  the  near-jet  region 
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Figure  22:  Instantaneous  planar  contours  of  the  streamwise  component  of  vorticity  for  the  steady  injection 
case. 
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Figure  23:  Instantaneous  planar  contours  of  the  streamwise  component  of  vorticity  for  the  pulsed  injection 
case  at  t/tp  =  0.5. 
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Figure  24:  Instantaneous  planar  contours  of  the  streamwise  component  of  vorticity  for  the  pulsed  injection 
case  at  t/tp  =  1.0. 
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List  of  Symbols 

Cp  =  turbine  blade  time-mean  surface  pressure  coefficient,  2(p  —  pj /  pitR 

V  =  duty  cycle  expressed  as  a  decimal  fraction 

Ek„ ,  Eu  =  nondimensional  turbulent  kinetic  energy  wave  number  and  frequency  spectra 

/  =  dimensional  jet  forcing  frequency 

J  =  Jacobian  of  the  coordinate  transformation 

kz  =  nondimensional  spanwise  wave  number 

l,n  =  nondimensional  blade  tangential  and  normal  directions 

M  =  Mach  number 

p  =  nondimensional  static  pressure 

Pr  =  Prandtl  number,  0.73  for  air 

q  =  nondimensional  planar  velocity  magnitude,  \J  v?  +  v2  or  \J v1  +  w2  for  the  jet 

Re oo  =  reference  Reynolds  number,  Poo9ooc/ Poo 

Re  =  chord  inlet  Reynolds  number,  ppi^Reaol /q 

t  =  nondimensional  time  based  upon  the  reference  velocity,  itoo 

tp  =  nondimensional  jet  forcing  period 

T  =  nondimensional  static  temperature 

in  the  x,y,z  directions 
«1,  lt2,  U3  =  u,  v,  w 

U,  V,  W  =  contravariant  velocity  components 

x, y,  z  =  nondimensional  Cartesian  coordinates  in  the  streamwise,  vertical,  and 

spanwise  directions 

xi, x2,x3  =  x,y,z 

a  =  turbine  blade  flow  angle 

7  =  specific  heat  ratio,  1.4  for  air 

<5  =  nondimensional  boundary-layer  thickness 

Sij  =  Kronecker  delta  function 

SV2,  =  2nd-order  and  4th-order  finite-difference  operators  in  £,77,  £ 

^4,  SV4,  6(4 

A l,  A?r,  A z  =  nondimensional  mesh  step  sizes  in  blade  tangential,  normal,  and  spanwise  directions 
p,  =  nondimensional  molecular  viscosity  coefficient 

■>  £,z  ■,  =  metric  coefficients  of  the  coordinate  transformation 

Vt,Vx,Vy,Vz, 

Ct )  Cx  j  C y  1  Cz 

to  =  nondimensional  frequency 

Subscripts 

i,o  =  inflow  and  outflow  conditions 

j  =  corresponding  to  vortex  generator  jet 

max,  min  =  maximum  and  minimum  values 

1  =  “local  freestream  value” 

00  =  dimensional  reference  value 

Superscripts 

+  =  value  given  in  law-of-the-wall  units 

—  =  time-mean  quantity 

'  =  root-mean-square  fluctuating  component 
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